Direct Detection of the Close Companion of Polaris with the 

Hubble Space Telescope^ 



O 

m 



Nancy Remage Evans, ^ Gail H. Schaefer,^ Howard E. Bond,^ Giuseppe Bono,'^ Margarita 
Karovska,^ Edmund Nelan,^ Dimitar Sasselov,^ and Brian D. Mason^ 

'66 

§ : ABSTRACT 
(N 

^ ■ Polaris, the nearest and brightest classical Cepheid, is a single-lined spec- 

troscopic binary with an orbital period of 30 years. Using the High Resolution 
Channel of the Advanced Camera for Surveys onboard the Hubble Space Tele- 
scope (HST) at a wavelength of ~2255 A, we have directly detected the faint 
1^ • companion at a separation of 0'.'17. A second HST observation 1.04 yr later con- 

I ■ firms orbital motion in a retrograde direction. By combining our two measures 

2 ■ with the spectroscopic orbit of Kamper and an analysis of the Hipparcos and FK5 

^ . proper motions by Wielen et al., we find a mass for Polaris Aa of 4.5!t^ 4 ^& — the 

first purely dynamical mass determined for any Cepheid. For the faint companion 
Polaris Ab we find a dynamical mass of 1.26^QQyM0, consistent with an inferred 
^ . spectral type of F6 V and with the flux difference of 5.4 mag observed at 2255 A. 

O ■ The magnitude difference at the V band is estimated to be 7.2 mag. Continued 

^ . HST observations will significantly reduce the mass errors, which are presently 

\Q , still too large to provide critical constraints on the roles of convective overshoot, 

^ . mass loss, rotation, and opacities in the evolution of intermediate-mass stars. 

O ■ Our astrometry, combined with two centuries of archival measurements, also 

confirms that the well-known, more distant (18") visual companion, Polaris B, has 
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a nearly common proper motion with that of the Aa,Ab pair. This is consistent 
with orbital motion in a long-period bound system. The ultraviolet brightness of 
Polaris B is in accordance with its known F3 V spectral type if it has the same 
distance as Polaris Aa,Ab. 

Subject headings: astrometry — binaries: visual — Cepheids — stars: funda- 
mental parameters — stars: individual (Polaris) 

1. Introduction 

Cepheid variable stars arc of central importance in galactic and extragalactic astronomy. 
They are the primary standard candles for measuring extragalactic distances, and they 
provide critical tests of stellar-evolution theory. Surprisingly, however, until now there has 
not been a single Cepheid with a purely dynamical measurement of its mass. 

Polaris {a UMi) is the nearest and, at 2nd magnitude, the brightest classical Cepheid, 
albeit one with a small light amplitude in its 3.97-day pulsation period (Turner et al. 2005 
and references therein). The amphtude, which had been decreasing for several decades, now 
appears to have stabilized and may be increasing again (Bruntt et al. 2008 and references 
therein) . The Hipparcos parallax of Polaris indicates a luminosity consistent with pulsation 
in the first overtone (Feast & Catchpole 1997). 

Polaris is the brightest member of a triple system (see Kamper 1996 and references 
therein). The well-known visual companion, Polaris B, is an 8th-magnitude F3 V star at a 
separation of 18". The Cepheid itself is a member of a single-hned spectroscopic binary with 
a period of 30 yr. In this paper we report the first direct detection of the close companion, 
from which we derive the first entirely dynamical mass measurement for a Cepheid. 

Cepheid masses are a key parameter for testing stellar evolutionary calculations. Be- 
ginning in the 1960's, discrepancies were found in the sense that Cepheid masses derived 
from pulsation modeling were lower than those derived from evolutionary tracks. A revision 
in envelope opacities brought pulsation and evolutionary masses closer together, partially 
alleviating this "Cepheid mass problem." However, recent evolutionary and pulsation con- 
straints for Galactic (Bono et al. 2001b; Caputo ct al. 2005; Keller 2008; Natale, Marconi, 
& Bono 2008) and Magellanic (Bono, Castellani, & Marconi 2002; Keller & Wood 2006) 
Cepheids still imply a discrepancy in masses at the ~15-20% level. The luminosities of 
Cepheids depend on the mass of the helium-burning core; physical mechanisms affecting the 
helium core mass include mixing due to convective core overshoot during the main-sequence 
phase, mass loss, stellar rotation, and radiative opacity. A directly measured mass for Polaris 
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would thus provide an important constraint on this theoretical framework. 

2. Observations and Data Reduction 

With the intention of a direct detection of the close companion, we imaged Polaris 
with the Hubble Space Telescope (HST) and the High Resolution Channel (HRC; plate scale 
0'/026pixel-^) of the Advanced Camera for Surveys (ACS). We chose the ultraviolet F220W 
filter (effective wavelength ~2255 A) in order to minimize the contrast between Polaris and 
the close companion, which wc anticipated to be a main-scqucncc star slightly hotter than 
the Cepheid, and also to minimize the size of the point-spread function (PSF). 

Observations were obtained on 2005 August 2-3, and again on 2006 August 13. At the 
first epoch we obtained a series of 0.1 to 0.3 s exposures dithered across 200 pixels on the chip 

over the course of one HST orbit, with several exposures taken at each dither position. At 
the second epoch we used the same dither pattern, but divided the spacecraft orbit between 
a series of 0.3 s exposures on Polaris A and 20 s exposures on Polaris B. For the longer 
exposures, we placed Polaris B at the same chip location as Polaris A in the short exposures, 
so as to provide an accurate PSF for a single star at the same place in the field of view. 

Figure 1 shows the co-added images of Polaris A from 2005 (left panel) and 2006 (middle 
panel). The close companion (which we designate Polaris Ab) is detected at the lower left 
of the primary (at about a "7 o'clock" position). Because of the asymmetric PSF shape, we 
performed several checks to confirm that the apparent companion is not an artifact. The 
right-hand frame in Figure 1 shows Polaris B in the 2006 image, with the star shifted to the 
same field position as Polaris A, and with its image scaled to the same flux as Polaris A. 
This PSF shows no artifact at the location of the companion seen in the images of Polaris A. 

Figure 2 shows contour maps of the same three images. Again the faint companion is 
seen in the 2005 and 2006 images of Polaris A, and there is no PSF artifact at this location in 
the image of Polaris B. We also retrieved images from the HST archive of several standard 
stars observed with ACS/HRC in the F220W filter over an interval of four years. Examples 
of these observations are shown in Figure 3. Although the PSF structure does vary somewhat 
with time, due to small changes in telescope focus and other instrumental phenomena, none 
of these images show any artifact at the location of the Polaris Ab companion seen in Figures 
1 and 2. 

To measure the separation and position angle between Polaris Aa and the newly revealed 
close companion, we used the calibrated flat-fielded exposures provided by the HST reduction 
pipeline. At each dither location, we median-filtered the repeated observations to remove 
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cosmic rays. We then extracted subarrays from the images centered on Polaris Aa with a 
size of (y.'85 x &!85. We used the observations of Polaris B from 2006 as a reference PSF to 
construct models of the close pair (Aa,Ab) by searching through a grid of separations and 
flux ratios. The IDL interpolate procedure was used to shift the PSF by sub-pixel intervals, 
using cubic convolution interpolation. The background was least-squares fitted with a tilted 
plane. 

The separations, position angles (PAs), and magnitude differences, determined through 
minimization between the models and the observations, are given in Table 1. The un- 
certainties were determined by analyzing multiple images individually and computing the 
standard deviation. We applied the filter-dependent geometric distortion correction of An- 
derson & King (2004) to convert the pixel values to a separation in arcseconds. To define the 
orientation of the detector y-axis on the sky, and thus determine the PA of the binary relative 
to the pole for the equinox J2000.0, we used the HST image-header keyword PA_APER. 

We also measured the separation and PA of the wide companion, Polaris B, relative 
to Aa, and the results are presented in Table 2. The good agreement with the historical 
measurements of the PA of Polaris B relative to A (see §4) indicates that we are properly 
defining the direction of north in spite of the extreme northerly declination. (Since the 
historical double-star convention is to give the PA for the equinox of the date of observation, 
we computed the precession corrections and give the adjusted PAs in the footnotes to Tables 1 
and 2, for the convenience of archivists.) 

3. Orbital Solutions 

We stress that the orbital analyses discussed below are only prehminary fits to data 
with a very limited sample (only two points) of separations and PAs. We followed three 
different approaches to determining the orbital parameters, in order to illustrate the scope 
of the available data. 

Kamper (1996) rederived the single-lined spectroscopic orbit of Polaris Aa using im- 
proved radial-velocity data, and a careful removal of the velocity signal due to the Cepheid 
pulsation. His solution provides the period, time of periastron passage, eccentricity, angle 
between the node and periastron, and the radial-velocity semi-amplitude of the primary 
(denoted P, T, e, a;, and Kas., respectively). By comparing the Hipparcos proper motion of 
Polaris Aa (which, over the duration of the Hipparcos mission, is nearly instantaneous in the 
context of the 30-year orbit) with the ground-based long-term average proper motion from 
the FK5 (which is essentially the motion of the center of mass) , Wielen et al. (2000) deter- 
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mined the inclination and the PA of the hne of nodes (i and Q). Their analysis, however, 
allows for retrograde and prograde orbital solutions (the two orbits being tangential at the 
Hipparcos epoch), with two different values of i and Q. The orbital parameters based on the 
Kamper (1996) and the two Wielen et al. (2000) solutions are presented in Table 3. 

The HST detection of the close companion Ab and its orbital motion at two epochs 

establishes a retrograde sense for the orbit (thus confirming the strong preference stated 
by Wiclcn et al. for their retrograde solution). Additionally, it provides constraints on the 
remaining unknown parameter of the orbit, the semi-major axis a. A combination of the 
spectroscopic mass function, 

/(M) = (MAbSini)V(MAa + MAb)' = 3.784 x 10-'kI^P{1 - , 



(1) 



, and the masses 

In §3.1-3.3 we describe the orbital fits that we computed based on a synthesis of the 
spectroscopic, astrometric, and HST data. In these sections, we present three successively 
more comprehensive orbital solutions. Table 4 summarizes the orbital parameters determined 
from each of these fits; the individual columns are described in more detail in §3.1-3.3. 

3.1. Semi-major Axis 

As a first approximation to an orbital solution, we fixed the spectroscopic and astro- 
metric parameters {P,T,e,uj,i,Q) to be those determined by Kamper (1996) and Wielen 
et al. (2000), and solved only for the semi- major axis a based on the two HST separation 
measurements. The orbital parameters for this solution are listed in the second column of 
Table 4. Figure 4 compares this retrograde orbit fit with the HST measurements, and shows 
extremely poor agreement for the PAs. However, the relatively large uncertainties in i and Q 
determined by Wielen et al. (2000) provide considerable flexibihty for adjusting the orbital 
parameters in order to improve the fit quality. 



with the total mass from Kepler s Third Law, 

Mtot = MAa + MAb = , 

then yields the masses of the binary components, 

M^. = 0.03357 ^^ty^ . 

MAa = Mtot - MAb , 

where a and the parallax tt are in arcseconds, P is in years, /^Aa is in km s 
are in Mq. 
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3.2. Best Fit to the HST Measurements 

To get a better fit to the HST data, we solved for 1], and a based on tlie two 
separation and PA measurements from tlie HST observations, wliile liolding tlie relatively 
well-determined spectroscopic parameters (P, T, e, and uj) fixed. We computed the orbit fit 
through a standard Newton-Raphson method in space, and present the results in the 
third column of Table 4. 

To explore the range of orbital parameters that fit the HST data, we performed a 
Monte Carlo search by selecting values of i, Vt, and a at random. We searched for 10,000 
solutions within the 3(j confidence interval, corresponding to a difference of Ax^ = 9 from the 
minimum value. Figure 5 shows cross-cuts through the surfaces for the three derived 
parameters. Using the recently revised Hipparcos parallax of 7.72±0.12 mas (van Leeuwen et 
al. 2007), we computed the total mass of the binary through Kepler's Third Law for all of the 
solutions found in the Monte Carlo search. In the last panel of Figure 5 we show a plot of the 
total mass versus inclination. Because a visual orbit is insensitive to the individual masses 
of the components, when combining the total mass with the spectroscopic mass function, 
there exist values of the inclination that produce negative masses for Mab. or MAb- In the 
remaining analysis we removed these negative-mass solutions from our sample of possible 
orbits. Essentially, this rejects all orbital solutions with i > 168°. The la uncertainties 
listed in the third column of Table 4 are determined from the Ax^ = 1 confidence interval 
of the modified distribution. The values of i and Q agree with the retrograde parameters 
computed by Wielen et al. (2000) at the l.So" and 2.2o" levels respectively. Figure 6 shows 
three examples illustrating how the orbit fit varies within a Icr confidence interval. 

We note that the HST-only solution yields a secondary mass of = 1-8 Mq (Table 4, 
col. 3), corresponding approximately to an A5 V star. The discussion below in §5.1, as well 
as the lack of a detection of the companion in UV spectra obtained by Evans (1988) with 
the International Ultraviolet Explorer {lUE ), makes it highly improbable that Ab could be 
this hot. With only two HST measurements sampling the orbit thus far, we do not yet have 
a good constraint on the curvature, and hence acceleration, of the orbit. In turn, this limits 
how well we can determine the total system mass, and it contributes to the large errors 
quoted in Table 4. 

3.3. Joint Fit to HST and Proper-Motion Measurements 

Incorporation of the Hipparcos proper-motion measurements into the orbital fit extends 
the time coverage of the measurements to ~15 years. This represents a significant fraction of 
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the orbital period and is therefore likely to improve the rehability of the results. Following 
the technique described in Wielen et al. (2000), we performed a simultaneous fit to the 
proper-motion data and the HST measurements. 

As Wielen et al. point out, the FK5 proper motion is averaged over several cycles of the 
orbital period and therefore reflects the center-of-mass motion of Polaris Aa,Ab. Because of 
the shorter time span of the Hipparcos mission, the Hipparcos proper motion more nearly 
represents an instantaneous measurement of the combined proper motion of the center of 
mass of the Aa,Ab pair and the orbital motion of the photocentcr about the center of mass 
at the epoch of the observations (~1991.3). The difference between the FK5 and Hipparcos 
proper motions thus gives the offset caused by the orbital motion. 

In computing the joint fit, we held the spectroscopic parameters {P,T,e,uj, and i^Aa) 
fixed at the Kamper (1996) values and solved for i, Q, and a, again using a Newton- Raphson 
technique in space. The input data were the relative positions of Aa and Ab at the 
two HST epochs, and the difference between the Hipparcos and FK5 proper motions. To 
incorporate the proper-motion data into the orbit fit, we had to compute the time-dependent 
offset of the photocenter relative to the center of mass predicted by the orbital parameters 
during the time of the Hipparcos observations. To compute these offsets, we converted the 
semi-major axis of Polaris Aa determined from the single-lined spectroscopic orbit to the 
semi-major axis of the photocentcr by using the mass ratio computed from the full set of 
orbital parameters and a magnitude difference between Aa and Ab of AV" = 7.2 (sec §5.1). 
This conversion is specified by cqns. (8) and (9) in Wielen et al. (2000). The expected 
difference between the instantaneous and the mean proper motion (A/igt) at the central 
epoch of the Hipparcos observations is then computed from eqn. (18) of Wielen et al. 

Table 5 shows the values of the proper motions used by Wielen et al. (2000). The first 
line shows the proper motion of Polaris given in the FK5 catalog (/UpKs)- In the second 
line, a systematic correction is applied to convert the proper motion from the FK5 reference 
system to the Hipparcos /ICRS system (Wielen et al. 2000). The proper motion measured by 
Hipparcos {/ihip) is given in the third line. The difference in the measured proper motions 
A/i — /ihip — A*FK5 is given in the fourth line. 

The last line in Table 5 shows the best-fit difference between the instantaneous and 
mean proper motion (A/xst) calculated from our simultaneous fit to the relative separation 
and PA measurements of Polaris Aa,Ab and the proper-motion data. Figure 7 shows a 
graphical representation of the Hipparcos and FK5 proper motions, and compares our best- 
fit value with the measured value for the difference between the instantaneous and mean 
proper motions. 
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As Figure 7 illustrates, the best-fit value of Afi in right ascension agrees well (O.Tcr) 
with the measured value, but the agreement is poorer (3.2cr) in declination. This discrep- 
ancy probably arises from our constraining the spectroscopic parameters to be exactly those 
derived by Kamper (1996), thus forcing the proper motions and HST measurements to ab- 
sorb the errors. We found that by allowing some variation in the spectroscopic parameters, 
we could substantially improve the fit to the proper-motion measurements. Once we have 
sampled enough of the visual orbit to better constrain i, Q, and a, the optimal orbital solution 
should be found by doing a simultaneous fit to the radial- velocity, HST, and proper- motion 
data. Unfortunately, however, Kamper only tabulated the radial velocities before removal 
of the pulsational variation, so a re-computation of the pulsation corrections would have to 
be carried out — a task well beyond the scope of this paper, and one that should await the 
availability of more HST observations. 

The orbital parameters and derived masses from our final combined fit are given in the 
last column of Table 4, which also contains our final best estimates for the dynamical masses 
of both stars. The best-fit orbit of Polaris Ab relative to Aa is plotted in Figure 8. The 
shaded-gray portion of the orbit marks the location of the companion during the interval 
of the Hipparcos observations in the early 1990 's, and it should be noted that the direction 
of motion at that time was, of course, 180° different from the direction of the differential 
motion of Polaris Aa (A/igt) shown as a green arrow in Figure 7. 



4. Astrometry of Polaris B 

Visual measurements of the PA and separation of the wide companion Polaris B relative 
to Polaris A extend back to the early 19th Century, with a few photographic observations 
being available from the 20th Century. We have compiled these measurements from Kamper 
(1996) and the Washington Double Star Catalog (Mason et al. 2001)0. 

In Figure 9 we plot the PA measurements of Polaris B, which have been precessed to 
the J2000.0 equinox. Due to its slow relative motion and large magnitude difference, Polaris 
A,B was generally not included in the major double-star observing programs, especially in 
the latter half of the 20th Century. This had the consequence that most of the measures that 
were obtained tended to be made by less-experienced observers, often using older techniques 
and/or smaller telescopes; this may explain the surprisingly large scatter in the late 20th 
Century. A linear least-squares fit to the data yields a rate of change in the PA of — 0?00035± 
0?00094 yr~^, consistent with no detectable change in PA for the past two centuries. (The 
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earlier work of Kamper had given a marginal detection of +0?0086 ± 0?0076 yr 

The separation measurements for Polaris B are plotted in Figure 10. There has been 
a slow downward trend in the separation, with a least-squares fit giving a rate of —1.67 ± 
0.19 mas yr~^. Kamper (1996), from the ground-based measures only, had found —1.7 ± 
0.6 mas yr~^. Since the absolute proper motion of Polaris A is ~46 mas yr~^ (see Table 5), 
the absolute motions of A and B agree to within about 4%. At the distance of Polaris 
given by the Hipparcos parallax, the difference in tangential velocities between A and B is 
1.03 ± 0.12 kms-^ 

In computing these least-squares fits, we weighted the observations by estimates of their 
measurement errors. For the HST observations, we applied the uncertainties quoted in 
Table 2. For the historical measurements, we divided the data into four groups spanning 
approximately 50 years each. We assumed measurement uncertainties equal to the standard 
deviation of the values measured in each of these four groups. 

The slowly diminishing separation of A and B (at constant PA) is not inconsistent with 
orbital motion in a physically bound pair — which is also supported by the close similarities 
of the radial velocities of A and B (Kamper 1996; Usenko Sz Klochkova 2007). To predict 
an order-of-magnitude rate of change in the separation, we assumed a circular orbit with a 
total system mass of M^a + + Mb = 7.15 Mq (based on the Aa-I-Ab mass of 5.8 
in the last column of Table 4 and a mass for Polaris B of 1.35 — see below). Adopting 
the revised Hipparcos parallax of 7.72 mas, and assuming an edge-on orbit with a period of 
~100,000 yr, we find a semimajor axis of a ^ 32" (or 0.02 pc). At the orbital phase implied 
by the observed separation of 18'.'2, the relative motion would then be —1.65 mas yr~^, close 
to the observed value. 



5. Astrophysical Properties of the Companions of Polaris 

5.1. Polaris Ab 

Our observed ultraviolet magnitude difference between Polaris Aa and Ab may be used 
to infer the spectral type, and hence the mass, of the newly resolved close companion. 

We downloaded UV spectra from the lUE data archiv^ for three F-type dwarfs having 



^The lUE data were obtained from the Multimission Archive at the Space Telescope Science Institute 
(MAST). Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant 
NAG5-7584 and by other grants and contracts. 
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accurate spectral types and parallaxes, as well as for Polaris itself. The F stars selected 
were 78 UMa (HR 4931, HD 113139; F2 V), HD 27524 (F5 V), and HD 27808 (F8 V). The 
latter two stars (Hyades members), as well as Polaris itself, were taken to be iinrcddcncd, 
while the spectrum of 78 UMa (a member of the Ursa Major group) was dcrcddcncd by 
E[B — V) = 0.01 mag. We then scaled the flux distributions for the three stars to the 
distance of Polaris, using the respective Hipparcos parallaxes. 

In Table 6, we list these comparison stars, their spectral types, masses implied by the 
spectral types, their Hipparcos parallaxes, absolute magnitudes based on the parallaxes, 
and flnally the predicted flux ratios relative to Polaris in the ACS/HRC F220W band. The 
adopted relationship between spectral types and masses is that of Harmanec (1988). The flux 
ratios were calculated by convolving the F220W system-throughput function (see Chiaberge 
& Sirianni 2007) with the scaled lUE spectra, and then ratioing with respect to Polaris. In 
Figure 11 we show the lUE spectrum of Polaris and the scaled spectra of the three F dwarfs. 

As listed in Table 1, the observed magnitude difference in the ACS/HRC F220W filter 
between Polaris Aa and Ab is 5.39 ± 0.08 mag, or a fiux ratio of 0.0070 ± 0.0006. A small 
correction to this ratio is needed because of the small (~10-15%) contribution to the signal 
from the red leak in the F220W filter, Polaris Aa being slightly redder than the companion. 

The red-leak contributions have been tabulated as a function of spectral type by Chiaberge 
& Sirianni (2007), leading to a corrected in-band flux ratio of 0.0074 ± 0.0006. This value 
is entered in the fourth row of Table 6, and is marked with a horizontal line in Figure 11. 
Interpolation in the last column of Table 6 then leads to an inferred spectral type of about 
F6 V, an absolute visual magnitude of My ~ +3.6, and an expected mass of 1.3 M0. 

The apparent y-band magnitude of Polaris Ab, for a distance modulus (m — M)o = 5.56, 
is inferred to be about 9.2, or some 7.2 mag fainter at V than Polaris Aa. This illustrates 
the advantage of observing the Polaris system in the ultraviolet, which lessens the contrast 
by nearly 2 mag. 

5.2. Polaris B 

As listed in Tabic 2, we also measured the F220W flux difference between Polaris Aa 
and B as 4.49 ± 0.04 mag, or a flux ratio of 0.0160 ± 0.0006. Correction for red leak, 
as described above, changes the flux ratio to 0.0178 ± 0.0006, entered in the second row of 
Table 6, and also marked with a horizontal line in Figure 11. Interpolating again in the table, 
we see that this ratio corresponds to a star intermediate between types F3 V and F4 V. In 
the case of the well-resolved Polaris B the optical spectral type has been determined from 
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the ground. Our result is in gratifying agreement with the spectral type of F3 V found by 
Turner (1977) and Usenko & Klochkova (2007), who also cite earlier spectral classifications 
of Polaris B by experts such as Bidelman. This finding not only validates our photometric 
analysis of Polaris Ab above, but again supports the physical association of Polaris A and B. 
Based on the relationship between spectral type and mass in Table 6, we infer the mass of 
Polaris B to be near 1.35 M©. 

Using the same method as for Polaris Ab, we can use the UV flux ratio to infer the 

V magnitude of Polaris B to be 8.7. The visual magnitude of Polaris B can be measured 
from the ground, but is made difficult by scattered light from Polaris A. Kamper (1996) used 
CCD imaging to determine a magnitude difference with respect to A of AV = 6.61, implying 

V — 8.59 (in good agreement with earlier photoelectric measurements oi V — 8.5 and 8.60 
by Fernie 1966 and McNamara 1969, respectively; Fernie included an approximate correction 
for scattered hght, and McNamara states that he observed only on excellent nights). In more 
recent work, to be reported separately, we have been carrying out astrometry of Polaris B 
with the Fine Guidance Sensors (FGS) onboard HST . As a byproduct, these observations 
yield an accurate V magnitude of 8.65 ±0.02. Thus our indirectly inferred V magnitude for 
Polaris B of 8.7 agrees very well with the ground- and HST-h&sed observations. 



6. Dynamical Masses 
6.1. Polaris Ab 

The final column in Table 4 lists the dynamical masses of both components of the 
close pair Aa,Ab obtained from our final orbital solution, as described in §3.3. For Ab, the 
dynamical mass is 1.26^o;o7M0. This is in remarkably good agreement with the 1.3 
inferred indirectly from the UV flux difference (§5.1), and is an indicator of the vahdity of 
our orbital solution. 



6.2. Theoretical Implications of the Cepheid's Dynamical Mass 

The dynamical mass of the Cepheid Polaris Aa from our flnal orbital solution, as listed 
in the last column in Table 4, is 4.5^^ 4 M©. 

We compare this result flrst with theoretical "evolutionary" masses, M^. The input data 
are the intensity-averaged mean apparent magnitudes (my = 1.98, uib = 2.58, from Fernie 
et al. 1995 and assumed to be unreddened), the revised Hipparcos distance of 129.5 ± 2.0 pc 
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(van Leeuwen et al. 2007), and a solar metal abundance (Luck &; Bond 1986; Usenko et 
al. 2005). We adopt the mass-period-luminosity (MPL) relation for He-burning fundamen- 
tal pulsators provided by Caputo et al. (2005, their Table 4). Before using this relation, 
we fundamentalized the first-overtone (FO) pulsation period of Polaris with the relation- 
ship logPp = logPpo + 0.13. By assuming Cepheid luminosities predicted by "canonical" 
evolutionary models that neglect convective-core overshooting, we find Me = 6.1 ± 0.4 M©. 
On the other hand, if we assume luminosities predicted by "noncanonical" evolutionary 
models that account for mild convective-core overshooting, given by L/Lcan — 1-3, we find 
Me = 5.6 ± 0.4 Mq. Use of the mass-color-luminosity (MCL) relation (Caputo et al. 2005, 
Table 5) yields very similar evolutionary masses. 

The HR diagrams in Figure 12 show a direct comparison between theoretical evolution- 
ary tracks and observations in the My, {B — V)o plane. In both panels of Figure 12 we plot 
the location of Polaris with an open triangle enclosing a small error bar. The top panel 
shows canonical evolutionary tracks at solar chemical composition; Me ~ 6M0 provides a 
good fit to the position of Polaris. The bottom panel shows noncanonical tracks, suggesting 
Me ~ 5.5 M0, except that the tip of the blue loop is not quite as hot as Polaris. However, 
the blueward extension of the loops is affected by chemical composition and by physical and 
numerical assumptions (Stothers & Chin 1991; Chiosi, Bertelli, & Bressan 1992; Bono et al. 
2000; Meynet & Maeder 2000; Xu & Li 2004). 

To compare our result with "pulsation" masses, Mp, we used the mass-dependent period- 
luminosity- color (PLC) relation of Caputo et al. (2005, their Table 2). We again funda- 
mentalized the pulsation period of Polaris, and using its intensity-averaged value of My 
find Mp(PLC) = 5.1 ± 0.4 M©. The pulsation mass of Polaris can also be estimated using 
the predicted period-mass-radius (PMR) relation for first-overtone Cepheids of Bono et al. 
(2001a), along with the radius of Polaris, it! = 46 ± 3 Rq, measured interferometrically by 
Nordgren et al. (2000). This gives Mp(PMR) = 4.9 ± 0.4 Mq. 

The lower pulsation masses, taken at face value, are thus in better agreement with the 
nominal dynamical mass than are the higher evolutionary masses. However, the current 
1(7 range of the measured dynamical mass, 3.1-6.7 M©, encompasses the entire range of 
theoretical masses. Thus, our discussion serves mainly to emphasize the crucial importance 
of reducing the error bars through continued HST high-resolution imaging of the Polaris 
Aa,Ab system. In addition, our companion HST FGS astrometric program will provide an 
improved trigonometric parallax. Doubtless there will also be future improvements in the 
spectroscopic orbit (e.g.. Turner et al. 2006; Bruntt et al. 2008). Simulations suggest that we 
can reduce the uncertainty on the Cepheid's dynamical mass to below ±0.6 Mq. This would 
provide a major constraint on the evolution of intermediate-mass stars and the physics of 
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Cepheid pulsation. 

6.3. Issues in the Evolution of Intermediate-Mass Stars 

As a further elaboration of the importance of accurate dynamical masses for Cepheid 
variables, we summarize the major open questions in the calculation of evolutionary tracks 
of intermediate-mass stars. A more complete discussion is provided by Bono, Caputo, & 
Castellani (2006). 

The luminosity of an evolved intermediate-mass star is related to the mass of the He- 
burning core. The physical mechanisms affecting the core mass include: 

1 . "Extra- mixing" of hydrogen into the core through convective core overshooting during 
the central hydrogen-burning phases; 

2. Mass loss, leading to a lower total stellar mass at the same luminosity; 

3. Rotation: the shear layer located at the interface between the convective and radiative 
regions enhances internal mixing, producing a larger He core mass; 

4. Radiative opacity: an increase in stellar opacity causes an increase in the central 
temperature, enhanced efficiency of central H-burning, and a higher core mass. 

Here we briefly discuss a few recent results that bear on convective overshoot and mass 

loss. 

The discussion of the Polaris mass in the previous subsection showed that the inclusion 
of noncanonical overshoot gave a better agreement with our preliminary dynamical mass. 
This is borne out by a mass measurement for the longer-period Cepheid S Muscae, based 
on HST Goddard High Resolution Spectrograph radial velocities of its hot companion, and 
an assumed companion mass based on its FUSE spectrum (Evans et al. 2006). The imphed 
mass of S Mus clearly favors mild convective overshoot. 

For mass loss, we note that the evolutionary calculations discussed in the previous 
subsection included semi-empirical mass-loss rates (Reimers 1975; Nieuwenhuijzen & de Jager 
1990). These rates are insufficient to resolve the discrepancy between evolutionary and 
pulsational masses. However, a variety of mostly recent observational information suggests 
that mass loss from Cepheids may be significant. At least two Cepheids, SU Cas and RS Pup, 
are associated with optical reflection nebulae (see Kervella et al. 2008 and references therein) 
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that may represent mass ejection from the Cepheids. Moreover, a large circumstellar envelope 
around the Cepheid ^ Car has been detected recently by Kervella et al. (2006), using mid- 
infrared data collected with the MIDI instrument on the VLTl. Circumstellar material has 
also been detected around 5 Cephei and Polaris itself by Merand et al. (2006). Recent Spitzer 
observations of a sample of Cepheids (Evans et al. 2007) have likewise revealed an infrared 
excess in the direction of 6 Cephei. 

On the theoretical side, a recent investigation (Neilson & Lester 2008) indicated that 
the coupling between radiative line driving (Castor et al. 1975) and the momentum input of 
both radial pulsation and shocks can provide mass-loss rates for Galactic Cepheids ranging 
from 10~^° to lO~^M0yr~^. This finding, together with typical evolutionary lifetimes (e.g.. 
Bono et al. 2000, Table 7), indicates that classical Cepheids may in fact be capable of losing 
the 10-20% of their mass that would be needed to resolve the discrepancy. 

7. Summciry 

The results of this study are as follows: 

1. We have used UV imaging with the ACS/HRC onboard the HST to make the first 
direct detection of the close companion of the classical Cepheid Polaris. 

2. We confirm orbital motion in a retrograde sense, based on two observations a year 
apart. 

3. By combining our HST measurements with the single-lined spectroscopic orbit (Kam- 
per 1996) and the FK5 and Hipparcos proper motions (Wielen et al. 2000), we derive a 
dynamical mass for the Cepheid Polaris Aa of 4.5^^ 4 — the first purely dynamical 
mass for any Cepheid. 

4. The dynamical mass is smaller than values estimated either from pulsational properties 
or evolutionary tracks, but the error bars are still large enough that the discrepancies 
have not achieved statistical significance. 

5. The close companion Polaris Ab has a dynamical mass of 1.26^^1^ Mq. This is consis- 
tent with a spectral type of about F6 V, inferred from the UV brightness of Ab. 

6. The more distant and well-known companion Polaris B has a UV flux consistent with 
its known spectral type of F3 V, lying at the same distance as the Cepheid. The 
proper motion of Polaris B is shown to be very similar to that of Aa,Ab, consistent 
with motion in a wide but bound orbit around the close pair. 
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7. Continued HST imaging, including two more observations that have been approved for 
our own program, will decrease the errors on the dynamical mass of Polaris, allowing 
a critical test of stellar-evolution theory and the influence of such effects as convective 
overshoot, mass loss, rotation, and opacities. 
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Polaris AaAb - 2005 




Polaris AaAb - 2006 




Polaris B --2006 




Fig. 1.— Coadded ACS HRC images of Polaris Aa,Ab taken with the F220W filter on 
2005 Aug 2 (left) and 2006 Aug 13 (middle). The close companion Ab is detected at the 
lower left of the primary (at about a "7 o'clock" position). The images are 0'.'85 x 0'.'85 and 
the directions of N and E are indicated. The right-hand panel shows a coadded image of 
Polaris B from longer exposures taken during the 2006 observations, and scaled to the flux 
level of the Polaris Aa,Ab images. There is no artifact in the Polaris B PSF at the location 
of Ab. 
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Fig. 2. — Contour maps of the coadded images shown in Figure 1. The outermost contour 
intervals in each panel range from 0.001 to 0.010 of the peak flux in steps of 0.001, and 
thereafter are at levels of 0.02, 0.04, 0.08, 0.16, 0.32, and 0.64 of the peak flux. The contours 
again demonstrate the absence of any artifact at the location of the Ab companion. 
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Fig. 3. — Archival images of flux-calibration standard stars observed with the ACS/HRC in 
the F220W filter. The star names and dates of observation are hsted in each panel. There 
is no artifact at the location of the Polaris Ab companion. 
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Fig. 4. — Comparison of the retrograde orbital solution of Wielen et al. (2000) (solid ellipse) 
with our HST measurements of Polaris Ab (filled circles). The arrow indicates the direction 
of motion. The orbit was calculated by fixing the spectroscopic and astrometric orbital pa- 
rameters (see Table 3) and solving only for the semi-major axis based on the HST separation 
measurements. We found a — 0'/131. The dotted lines connect the observed positions to 
x's marking the predicted positions, and show the inadequacy of this simple solution at 
predicting the position angles. The size of the HST error bar is indicated by the cross at 
the top of the plot. 
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Fig. 5. — Crosscuts through the surfaces derived from the fit to the HST measure- 
ments alone while holding the spectroscopic parameters fixed. These results were obtained 
by performing a Monte-Carlo search for orbital solutions within the 3a confidence interval 
(Ax^ = 9). The total mass was derived for each of the 10,000 orbits found within this inter- 
val using a parallax of 7.72 mas. The dotted lines indicate the critical value of the inclination 
(z = 168°), above which negative values of a component mass are produced. The color codes 
in the electronic version of this figure correspond to the Icr (red), 2a (blue), and 3a (green) 
confidence intervals. 



-22- 



0.1 



0.0 - 



Q_ 
CD 



0.1 



-0.2 



T — I 1 1 1 1 1 — I r 



T 1 1 — I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 

+ HST error bor 




2005.59 
2006.62 



i=164°, a=0.115", Q = 7.5°, Mtot=3.6M0 
i=155°, a=0.116", = 9.0°, Mtot=3.9M0 
1=139°, a=0.125", Q=14.°, Mt^t=4.9Mo 



J I I I I I I I L 



J I I I I I I I L. 



J I I I I I I I L 



0.1 



0.0 -0.1 

E-W sep (") 



■0.2 



Fig. 6. — Three examples of orbits resulting from variations of the orbital parameters within 
the 1(7 confidence intervals based on the fit to the HST measurements alone while holding 
the spectroscopic parameters fixed. The arrows indicate the direction of motion. All three 
orbits fit the HST measurements (filled circles) within the la error bar shown at the top of 
the figure, but they imply total system masses ranging from 3.6 to 4.9 M©. 
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Fig. 7. — Schematic comparison between the instantaneous and mean proper motions of 
Polaris A. The vectors indicate the magnitude and direction of the annual proper-motion 
measurements from Hipparcos (/xhip; red), the long-term average FK5 ground-based motion 
(/Wfk5; blue), and the vector difference between the two (A/i; black). Dotted red and blue 
lines indicate the la uncertainties in the Hipparcos and FK5 measurements. The green 
vector represents the best-fit difference between the proper motions (A/ifit) computed from 
our simultaneous orbit fit to the HST measurements and the proper-motion data while 
holding the spectroscopic parameters fixed (see §3.3). 
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Fig. 8. — Final best-fit orbit of Polaris Ab relative to Aa based on a simultaneous fit to 
the HST measurements (filled circles) and the proper- motion data while holding the spec- 
troscopic parameters fixed. The predicted orbital positions at the epochs of the HST mea- 
surements are marked by the x symbols (partially hidden by the observed points). The 
position of the companion during the time-frame of the Hipparcos mission is highlighted by 
the shaded gray line segment; its direction of motion is of course 180° different from the 
direction of A/Xfit, shown as a green arrow in Figure 7. 
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Fig. 9. — Historical measurements of the position angle of Polaris B relative to A, precessed 
to the equinox of J2000.0. Our HST measurements are marked by the two filled circles. 
The solid line shows a least-squares fit yielding a rate of change in the position angle of 
-0?00035 ± 0?00094 yr-^. 
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Fig. 10. — Separation measurements of Polaris B relative to A. Our HST measurements are 
marked by the two filled circles. The sohd line shows a hnear least-squares fit yielding a rate 
of change in the separation of —1.67 ± 0.19masyr~^. 
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Fig. 11. — lUE UV spectra of Polaris, and of three main-sequence F-type stars scaled to 
the distance of Polaris using the respective Hipparcos parallaxes. The two horizontal lines 
correspond to the mean flux levels in the F220W fllter of the ACS/HRC of Polaris B (top) 
and Ab (bottom). The lengths of the horizontal lines correspond to the FWHM of the filter. 
On the basis of its UV fiux level, Polaris B is inferred to have a spectral type near F3-F4 V, 
in good agreement with its ground-based classification at F3 V. Polaris Ab is inferred from 
its F220W fiux to have a spectral type near F6 V. 
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Fig. 12. — Evolutionary tracks (see Pietrinferni et al. 2006) in the My vs. {B — V)o plane, 
with the location of Polaris marked {open triangle enclosing error bar). Top panel: canonical 
evolutionary models neglecting convective core overshooting during H-burning phases. Bot- 
tom panel: noncanonical evolutionary models including convective core overshooting. The 
assumed Polaris distance and hehum (Y) and metal (Z) abundances (scaled solar) are in- 
dicated in the figure, and the legends on the left indicate the stellar masses for each track. 
These models include mass loss with a Reimers parameter rj — 0.4. 
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Table 1. Position and magnitude of Polaris Ab relative to Polaris Aa 



Besselian Date 


UT Date & Time 


Pi") 


PA (J2000) i")"" 


Am(F220W) 


2005.5880 


2005 Aug 2 23:45 


0.172 ± 0.002 


231.4 ± 0.7 


5.38 ± 0.09 


2006.6172 


2006 Aug 16 22:01 


0.170 ± 0.003 


226.4 ± 1.0 


5.40 ± 0.09 



^PAs for equinox of date are 232?8 and 228? 1. 



Table 2. Position and magnitude of Polaris B relative to Polaris Aa 



Besseliaii Date 


l^T Date t Time 


p n 


PA (J2000) (°)'"' 


Ar»(F220AV) 


2005.5880 


2005 Aug 2 23:45 


18.217 ± 0.003 


230.540 ± 0.009 


4.53 ± 0.04 


2006.6172 


2006 Aug 16 22:01 


18.214 ± 0.003 


230.520 ± 0.009 


4.45 ± 0.02 



^PAs for equinox of date are 231?980 and 232?216. 



Table 3. Previously determined orbital parameters 



Kamper (1996) 


Wielen et al. (2000) 
Prograde 


Wielen et al. (2000) 
Retrograde 


P (yr) 29.59 ± 0.02 
T 1987.66 ± 0.13 

e 0.608 ± 0.005 
cu (°) 303.01 ± 0.75 
XAa (km s-1) 3.72 ± 0.03 


i (°) 50.1 ± 4.8 
(°) 276.2 ± 9.5 


i (°) 130.2 ± 4.9 
(°) 167.1 ± 9.4 



^Values of fl quoted in Wielen et al. (2000) correspond to the astrometric orbit 
of Polaris Aa relative to the center of mass. 
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Table 4. Orbital Parameters^ and Masses^ 



Parameter Wielen et al. (2000) Fit HST Joint Fit to HST 
Retrograde^ only & Proper Motion 



in 


130.2 (F) 


15511^ 


128l^t 


Vt^ (°) 


347.1 (F) 


9.011^ 


19lf 


«(") 


0.131 ±0.04 


0.1161^°^ 


n 1 OO+0.015 
U.100_o oil 


Mtot (Mq) 


5.6 ± 1.0 


Q n+io 

■^■^-0.5 


5.8i?-2 


MAa (Mo) 


4.3 ± 1.3 


^■-'--2.1 




MAb {Mq) 


1.26 ±0.80 




1 26+°-^'' 



^Spectroscopic orbital parameters (P, T, e, a;) held fixed at Kamper 
(1996) values. 

''Assumes Hipparcos parallax of 7.72 ± 0.12 mas. 

^Orbital parameters marked (F) are fixed at the given values when 
computing the best-fit solution. 

'^Q has been rotated by 180° from the values quoted in Wielen et al. 
(2000) to correspond to the orbit of Polaris Ab relative to Polaris Aa. 



Table 5. Proper-motion Data for Polaris A 



Quantity 


System 


Ha cos 6 


fJ'S 






(mas yr~^) 


(mas yr~^) 




FK5 


+38.30 ± 0.23 


-15.20 ±0.35 


A*FK5 


HIP 


+41.50 ± 0.97 


-16.73 ±0.75 


A*HIP 


HIP 


+44.22 ± 0.47 


-11.74 ±0.55 


A/i = //hip — A*FK5 


HIP 


+2.72 ± 1.08 


+4.99 ± 0.93 


A/ifit 


HIP 


+4.59 ± 2.52 


+1.21 ± 0.74 
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Table 6. Spectral- Type Comparisons 



Star 


Spectral 


Mass 


Parallax 


My 


F220W Flux 




Type 


(Mo) 


(mas) 




Relative to Polaris'^ 


78 UMa 


F2 V 


1.41 


40.06 


+2.9 


0.0262 


Polaris B 










0.0178 


HD 27524 


F5 V 


1.33 


19.55 


+3.2 


0.0114 


Polaris Ab 










0.0074 


HD 27808 


F8 V 


1.22 


24.47 


+4.1 


0.0029 



'^Flux ratios for 78 UMa, HD 27524, and HD 27808 are predicted from 
their IJJE spectra, scaled to the distance of Polaris; ratios for Polaris B 
and Ab are those observed by us, corrected for red leak as described in 
the text. 



